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Abstract 

We study the large-distance behavior of the bulk order-parameter 
correlation function G(r) for T > T c within the lattice version of the y? 4 
theory including lattice effects. We also study the large- L behavior of 
the susceptibility \ f° r T > T c of the confined lattice system of linear 
size L with periodic boundary conditions. We find that the structure 
of the large- L behavior of x of the confined system is closely related to 
the structure of the large-distance behavior of G(r) of the bulk system. 
Explicit results are derived in the spherical (large-n) limit and in one- 
loop order for general dimensions d > 2. For the lattice model with 
cubic symmetry we find that finite-size scaling must be formulated in 
terms of the anisotropic bulk correlation length (exponential correla- 
tion length) that governs the exponential decay of G(r) for large r 
rather than in terms of the ordinary isotropic bulk correlation length 
£ defined via the second moment of G(r). We show that it is the ex- 
ponential bulk correlation length £i in the direction of the cubic axes 
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that determines the exponential finite-size scaling behavior of lattice 
systems in a rectangular geometry. This result modifies a recent in- 
terpretation concerning an apparent violation of finite-size scaling in 
terms of the second-moment correlation length £ ^ £1. Exact results 
for the one-dimensional Ising model illustrate our conclusions. Fur- 
thermore we show for general d > 2 that a description of finite-size 
effects for finite n in the entire region < L/£ < oo requires two dif- 
ferent perturbative approaches that are applicable either to the region 
< L/£ < 0(1) or O(l) < L/£ < oo, respectively. In particular we 
show that the exponential finite-size behavior for L/£ » 1 above T c 
is not captured by the standard perturbation approach that separates 
the homogeneous lowest mode from the inhomogeneous higher modes. 
Consequences for the theory of finite-size effects above four dimen- 
sions are discussed. We show that the two-variable finite-size scaling 
form predicts an exponential approach oc e~ L ^ x to the bulk critical 
behavior above T c whereas the reduction to a single-variable scaling 
form implies a power-law approach oc L~ d . 
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1 Introduction 



The fundamental length scale near a critical point is the bulk correlation 
length £ which is a measure of the range of correlations of the order-parameter 
fluctuations. In this paper we consider lattice systems with cubic symmetry 
at and above T c . We assume short-range interactions. The bulk order- 
parameter correlation function G(r) between two lattice points at relative 
separation r serves to define £ via the second moment of G(r) according to 

£ 2 = ^r 2 G(r)/^G(r) . (1.1) 

r r 

In terms of the Fourier transform G(k) this definition reads 

^GtOjA^k)]-^. (1 . 2) 

This is applicable to lattice systems with cubic symmetry whose correlation 
function G(k) has an isotropic small k behavior at 0(k 2 ). The "second 
moment correlation length" (|1.2p is widely used in field-theoretic calculations 
[jl], ^ 0] since £ is identical with the inverse mass m~ l = £ that enters the 
two-point vertex function r( 2 )(k) = G(k) -1 . 

The bulk correlation length ( |1.2|) plays a fundamental role also in the for- 
mulation of the finite-size scaling behavior of confined systems. Consider, 
for example, the susceptibility x{t, L) oc 6(0) of a ferromagnetic system for 
t = (T — T c )/T c > in a finite geometry with a characteristic size L and with 
periodic boundary conditions. It is believed that the relative deviation from 
bulk critical behavior x(t, oo) = A x t~' y has the asymptotic (large L, small t) 
scaling form @, ||, ||, 0, || below four dimensions 



Ax = L) = 9W0 (1.3) 



where g(x) is a universal function in the entire range < x < oo. This 
means that any dependence of A% on the lattice constant a is negligible in 
the asymptotic region £ ^> a, L ^> a for arbitrary ratios including the 



large- L behavior at fixed £ ^> a j|, |9| . In a recent paper we have shown [10 
that this statement is not valid in the regime L/£ ^> 1. Specifically, for cubic 
geometry of size L 3> £ ^> a, we have found 



A X = s(L/Oexp jr(a/0|| • (1.4) 

Here g(x) is indeed universal but the a dependent exponential factor with 
the nonuniversal function 



na/0 = y 4 (a/0 2 + O[(a/0 3 } (1.5) 



is non-negligible for L > 24 £ 3 /a 2 , i.e., for sufficiently large L close to the 
bulk limit above T c , even in the asymptotic region £ ^> a. Thus the finite-size 
scaling form ( |1.3|) in terms of the second-moment correlation length (|1.2|) is 
not valid for lattice systems, even below four dimensions. This conclusion is 
based on one-loop results as well as on exact results in the spherical limit of 



the if lattice model at finite lattice spacing [10]. We note that the possibility 



of a non-negligible dependence of finite-size effects on the lattice spacing a, 
even for d < 4, was already mentioned by Privman and Fisher p| . An as 



yet unexplained a-dependence that is not consistent with a finite-size scaling 



ansatz of the kind (|1.3|) was also found by Gelfand and Fisher [12] in the 



interfacial free energy of the two-dimensional Ising model for T < T c 



No intuitive reasoning was given in our recent work [|T0| as to what might be 
the physical origin for this unexpected failure of the finite-size scaling prop- 
erty. It is the purpose of the present paper to elucidate this unsatisfactory 
situation by further analyzing the role played by the correlation length in 
both bulk and confined lattice systems. Instead of the second-moment cor- 
relation length (|1.2p we consider the "exponential correlation length" £ e that 
governs the large- distance behavior of the bulk system at fixed T above T c 

ISA, 



G{re) ~ B d r^ 2 exp(-r/£ e ) , (1.6) 

where e = r/r is the unit vector in the direction of r and where £ e is defined 
by 



C^-lim {[In G(re)]/r}. (1.7) 



We note that ( |1.6j ) is expected to hold even well above T c outside the asymp- 
totic critical region. For lattice systems, £ e is an anisotropic quantity. It is 

5 



expected, however, that the asymptotic ratio lim f ^ 0+ £ e /£ becomes isotropic 
and has a universal value fl4|, [15], [16], [17], Hi, in agreement with our results. 



In the present paper we shall call attention to the fact that the non-universal 
difference between £ and £ e is non-negligible even asymptotically close to T c 
if the scaling function has an exponential form. We shall show that lattice 
effects cause additive non-universal corrections to the asymptotic form of £ e 
[see ( |1.9| ) below]. Since £ e appears in the exponent of the scaling functions 
these additive corrections turn into non-negligible multiplicative overall fac- 
tors for the exponential scaling functions of both G(r) for r ^> £ e and x(t, L) 
for L > £ e . 

Because of the anisotropy of £ e for lattice systems, there exist infinitely many 
£ e in contrast to the unique isotropic quantity £. Since x(A L) oc G(0) oc 
G(r) involves all directions of r it is not clear a priori whether a certain 
average of £ e or £ e in a particular direction e enters the finite-size scaling 
form of x{t,L). Here we shall find that for cubic geometry the exponential 
finite-size effect is determined by the correlation length £i = £ (and not by 
£ ^ £i) where ei is the unit vector in the direction of one of the cubic axes. 
Specifically we derive an explicit relation between £i and £ in the spherical 
limit (Sect. 3) and in one-loop order of the (p A lattice model (Sect. 4). This 
relation reads at finite lattice spacing a 




(1.8) 
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It is remarkable that the non-universal a dependence of (|1.4j) is completely 
absorbed by the exponential correlation length £1 if the function £(£i) of ( |1.8| ) 
is substituted into the right-hand side of ( |1.4j ). In the present paper we shall 
show that for L 3> £ a 



^(L/e(6))e^|r(a/e(ei))^} =P(V6) , (1.10) 

hence universal finite-size scaling below four dimensions at finite lattice spacing 
is restored in the form 



A x = g(L/^) (1.11) 
with the universal function g(x). 

Thus, at one-loop order and in the spherical (large-n) limit, there is no 
violation of finite-size scaling at finite lattice spacing below four dimensions 
provided that the exponential correlation length £ x rather than £ is employed 
as the bulk reference length. We note that this result is not a general con- 
sequence of the renormalizability of the <£> 4 theory but is only an explicit 
computational result for cubic (and rectangular) geometry in one-loop order 



and in the spherical (large-n) limit. It is also valid for the one-dimensional 
Ising model (Sect. 6). At present it is an open question whether this result 
remains valid at finite n beyond one-loop order and whether it can be based 
on more general arguments. The renormalization-group (RG) arguments in 
the limit a — > || are not sufficient to establish the complete finite-size 
scaling form for lattice systems. Although these RG arguments for a — > 
lead to the same universal scaling function g(x) as our analysis at finite a, 
they do not identify the a dependent finite-size scaling variable x = 

We shall also show (Sect. 4) that two different perturbative treatments are 
necessary to describe the finite-size effects in the entire asymptotic region 
< L/£ < oo and that the previous finite-size perturbation approach below 



four dimensions [19L EQ, 21] does not capture the exponential structure of the 



scaling function g{x) oc e x for x ^> 1. 

The necessity of employing £i rather than £ is not restricted to d < 4 dimen- 
sions. In Sect. 5 we shall discuss the consequences of our results for lattice 
systems with d > 4 where £j rather than £ should be incorporated in the 
two- variable finite-size scaling form [fj2|, ^3f. The inhomogeneous modes are 
shown to yield Ax oc e _L ^ 1 for L ^> £ whereas the lowest-mode approxi- 
mation |T!| and the corresponding single- variable finite-size scaling form |2| 
imply Ax oc L~ d . 

Very recently an exponential size dependence has been confirmed by Stauffer 
p5| in Monte Carlo simulations for the magnetization of the Ising model in 
two, three and five dimensions . 
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In summary, even though we will confirm through ( |1.9j ) (within our approx- 
imations) that asymptotically close to the critical point the two correlation 
lengths £ and £1 are the same and isotropic, the scaling form g(x) of the 
leading finite-size effect near the bulk limit is sensitive to the choice of the 
correlation length because of the exponential dependence of g(x) oc e~ x on 
the correlation length. 
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2 Lattice effect on the large-distance behav- 
ior of the bulk order-parameter correlation 
function 

We consider a y? 4 lattice Hamiltonian for the variables ipi on the lattice points 
Xj of a simple-cubic lattice with the lattice constant a. The variables tfi have 
n components ip ia with a = 1, 2, n which vary in the range — oo < <p ia < 
oo. We assume the statistical weight oc e~ H with 



where are the dimensionless couplings of a short-rang interaction with 
cubic symmetry and where fee T w fc B T c is absorbed in r ,-u and J^. The 




(2.1) 



variables v?j have the dimension [a^ 2 d )/ 2 ] . We are interested in the large- 



distance behavior of the bulk correlation function 




(2.2) 



k 



above T c , normalized to the number of components n, where 




(2.3) 



j 
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with some fixed lattice point x . In ( ^.2[ ) J k stands for (2ir)~ d J d d k with a 
finite lattice-cutoff | k m \< 7r/a,m = 1,2, .. . , d. First we consider the limit 
n — ► oo at fixed won for d > 2 in which case we obtain |p7] 



G(k)- 1 = G(O)- 1 + J k 



(2.4) 



J k = -[J(0)- J(k)], 



(2.5) 



J(k) = (a/L) d ^ Jij e~ ik ^-*J . 



(2.6) 



G(0) is determined by an implicit equation |27| which, however, will not 
be needed in the following since G(0) can be expressed directly in terms of 
£ 2 . Using the second-moment definition for the bulk correlation length £ 
according to (|1.2j ) we have |27 



G(0) = Jo" 1 ^ 



(2.7) 



Jo = ^(S/ J L) d ^(VS 2 )(x l -x, 



(2.8) 
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For simplicity we consider a nearest-neighbor interaction J > which yields 



J(k) = 2 J ^2 cos{k m a), (2.9) 



m=l 



4J ' 



k = — - cos(A; m S)] = J k 2 + 0(k?kj) (2.10) 



a 



with Jq = 2J. 

In summary we need to calculate the large-distance behavior of 




G(x) = / e H (S/0 2 + 2 £[1 - cos(£; m a)] (2.11) 



m=l 



where we have chosen xo = and = x = (xi, X2, x,i) with Cartesian 
coordinates x m . Eq. (|2.11 ) is valid not only in the large-n limit but also in 



an ordinary perturbation calculation to one-loop order for general n. (For 
the latter renormalization-group treatment is carried out in Appendix 

B for d < 4.) A representation of G(x) in terms of Bessel functions of integer 
order v (see, e.g., 9.6.19 of Ref. |2£ 



7T 

I v {z) = i y d0 cos9 cos(v6) (2.12) 
o 
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can be given as 



G(X) = 2J 



ds e 



-(a/0 2 s 



e ikx exp 



-2s^2[l-cos(k m a))\ (2.13) 



m=l 



oo 

a 2 - d 



2J 



I ds e'^ s e' 2ds f[ I Vm {2s) (2.14) 
o m=1 



with the integers v m = x m /d . In general, G(x) is an anisotropic function 
whose exponential large-distance behavior ~ exp (— |x|/£ e ) leads to the defi- 
nition of an anisotropic correlation length £ e ^ £ in the direction of the unit 
vector e = x/|x| [I3|. An explicit demonstration of the anisotropy of £ e is 



given in Appendix A where the angular dependence of £ e is calculated for 
the case where e lies in the 2-dimensional x\ — X2 plane of the d-dimensional 
bulk system. For our present purpose it suffices to consider only the special 
case where e = ei = (1, 0, 0...) is the unit vector along one of the cubic axes. 
Then we have x = (x, 0, 0...) and kx = k x x . The corresponding correlation 
function is denoted by C(x) = G(x) which is obtained from ( |2.14j ) as 



C(x) = \ T j ds e~^ 2s e- 2ds [h{2s)} d - l I xra {2s) . (2.15) 
o 

This result is valid for arbitrary x/d and £/a and therefore does not yet have 
a scaling form. In Appendix A the large- |x| behavior of C(x) at arbitrary 
fixed £/o is derived. The result is 
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4J \2ir\x\J 

We see that in the large- \x\ limit the 
tial correlation length £i = £ ei in the 
than the second-moment correlation 
and £ for n — > oo reads 



(2.16) 

natural reference length is the exponen- 
direction of one of the cubic axes rather 
length £. The exact relation between £1 



T^-sinh (A 



(2.17) 



or 



(2.18) 



The difference between £ and ^ is a true lattice effect that disappears in 
the formal limit a — > 0. Eqs. Q2.17| ) and ( 2.18 ) are also valid for finite n in 
one-loop order above T c but in two-loop order and beyond we expect (small) 
corrections to ( |2.17| ) and ( |2.18 ) for finite n. 



Eqs.( |2.16] )-( |2.18| ) are valid for arbitrary £i/a even well above T c . In the 
asymptotic region £i 3> a, (|2.16| ) attains the scaling form 
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C{x) ~ (a/\x\) d - 2+ v $(1x1/6) 



(2.19) 



with the scaling function for |x|/£i 3> 1 



$(N/6) = A 



a 2 ~ d 



U(2n f- 1)/2 V6 



i(d-3)+r? 



exp (-|x|/fi 



(2.20) 



where rj = and A = 1 in the present case of the limit n — > oo. Eqs. 
(|2.19p and ( |2.20|) are also valid for general n in one- loop order for d < 4, see 
Appendix B. In this case we have a critical exponent r\ > and an amplitude 
i / 1 which we obtain from a RG treatment at finite a ||10|| , applied to the 
bare one-loop result ( p.llQ for general n, as described in Appendix B. Eqs. 
2.19| ) and ( |2.20 ) are also valid in one-loop order for general n and d > 4 



where rj = and A = 1, apart from O (mq) corrections. 



Close to T c where both £ and £i diverge, an expansion of ( 2.1 8|) yields 



6 = £ 



i + 



24 U 



~\ 2 

a 



+ 



(2.21) 



Thus, for n — > oo, £i and £ become identical sufficiently close to T c . (This 
is also valid for finite n in one-loop order above T c but in two-loop order 
and beyond the asymptotic value of the ratio for T — > T c is expected to 
become different from 1, see e.g. Refs. [jT4 , |i~5l , [16| , |I~7l , 0.) Therefore £i can 
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be replaced by £ in the prefactor of ( |2.20| ). We emphasize, however, that a 
replacement of £1 by £ is not possible in the exponential part of ( [2.16D and 



( |2.20| ), even arbitrarily close to T c . This is seen by substituting ( |2.21| ) into 
the exponential function of ( [2.16| ) and ( |2.20| ), 



e 



-|x|/{i = exp (| x |a 2 /24£ 3 ) e" N/? , (2.22) 



C(x) ~ (a/M) d-2+7? $ (\x\/£) exp (|x|a 2 /24£ 3 ) . (2.23) 



Now the additive correction in ( |2.21| ) has turned into an exponential non- 



universal prefactor in Q2.22 ) and (|2.23| ) that cannot be simply replaced by 1 
and that is by no means negligible for sufficiently large \x\ > 24£ 3 /a 2 , even 
in the asymptotic critical region £ 3> a. This is the crucial point of our 
argument. 

Thus, in order to have a universal (a independent) scaling form of C(x) 
for large \x\ ^> £ at fixed T above T c where C(x) has an exponential form, 
it is inescapable to employ £i rather than £ as the appropriate reference 
length. Correspondingly, for any T > T c , there exists an infinitely large 
region \x\ > 24 £ 3 /a 2 where the anisotropy of G(x) is no longer a negligible 
correction to the isotropic part. In the critical region a |x| £, on the 
other hand, where G(x) has a power-law behavior the nonuniversal part of the 
difference between £ and £ e can be considered as a negligible non-asymptotic 
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additive correction. The natural reference lengths in these two regions are £ 
and £ e , respectively, and a complete scaling description should embody a kind 
of crossover in the scaling variable from |x|/£ to x/£ e . The same situation 
will arise in the finite-size problem with respect to the L dependence of the 
susceptibility that is analyzed in the next Section. 

The analysis of this Section can be extended to the continuum version of the 
(f 4 theory. In a separate paper |^9[ we shall show that the results depend on 
the cutoff procedure. An (isotropic) exponential large- |x| behavior of G(x) 
is found for a smooth cutoff whereas a sharp cutoff implies a nonuniversal 
non-exponential |x|-dependence of G(x). In the case of a smooth cutoff it 
is also found that the exponential correlation length differs from the second- 
moment correlation length even though the continuum system is isotropic. 
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3 Lattice effect on the finite-size scaling be- 
havior for n —> oo 



We consider the lattice Hamiltonian (2.1) for a finite hypercubic geometry 
with volume V = L d and with periodic boundary conditions. We are inter- 
ested in the exact large-L behavior of the susceptibility x — G(0) above T c 
in the large-n limit at finite lattice spacing a. Specifically we wish to identify 
the reference length that governs the expected exponential L-dependence at 
fixed T > T c . The answer is not clear a priori since in the sum 



a d 



X=— y2<(pj(p > (3.1) 



n 

3 



there are contributions from < ifj (po > in all directions involving all anisotropic 
correlation lengths £ e discussed in the preceding Section. The finite-size ef- 
fect on x a t finite lattice constant a has already been calculated previously 



T0| , |22| , |27|| where it was expressed in terms of the second-moment correlation 
length £. Here we shall demonstrate that £i as calculated in the preceding 
Section, rather than £, is the appropriate reference length in the finite-size 
scaling structure. 



We start from the implicit equation for n — > oo at fixed uqti [27 



X 1 = ro 



+ 4u nL- d J2(^ + X' 1 ) 1 (3.2) 
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which can be rewritten as 



x 1 = r -r 0c + Au nD(x 1 ,L,a) 



4u n x 



-i 



4 4 + x 



-i 



-i 



(3.3) 



where r 0c = — 4-u n J k J k . The finite-size effect is contained in the function 



D(x~\L,a) =L- d J2(-J* 



X 



-i 



-i 



-i 



(3.4) 



d~s e ~~ s/x I L- d 



(3.5) 



The summations run over discrete k vectors with components kj = 2irmj/L, rrij 
0, ±1, ±2, = 1, 2, d, in the range —it /a < kj < n/a. Since is a pe- 
riodic function of each component kj the sum in ( |3.5| ) satisfies the Poisson 
identity g, || 



(3.6) 



where kn = £V kjUj. The sum runs over all integers rij,j = 1,2, 

in the range — oo < rij < oo whereas and J k have finite cutoffs ±7r/a. 

For the case of a nearest-neighbor coupling J > we have 
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-sJk 



n exp 



m=l 



4Js 



[1 - cos(/c m a)] 



(3.7) 



This leads to the representation in terms of the Bessel functions I v (z), ( |2 . 1 2| ) 



-S.Ju 



[a- 1 e- 2s I (2s)}' 



(3.8) 



-SJ k 



1 £ e ~ 2s ^( 2 *) 



(3.9) 



where s = 2sJ/a 2 and where u n = nL/a are integers. The resulting expres- 
sion for D reads 



D( X ~\L,a) 



d 2 ~ d 



2J 



ds e - a V(2Jx) e - 2ds 



/ (2 S ) + E 2 ^(2 S ) 



n=l 



- [I (2s)} d \ . (3.10) 



The large- L limit corresponds to large integers v n = nL/d. Since we consider 
this limit at fixed temperature above T c we may replace x(tj L) in the expo- 
nent of ( |3.10|) by the bulk value Xb — ^ 2 /2J. For the asymptotic behavior of 
I u (2s) we refer to Appendix A. The leading term for L ^> £ comes from the 
n — 1 contribution in (|3.10|) , 
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D( X -\L,~a) = {d ~a 2 - d /J) j dse-^ 2s e~ 2ds [Ws)]^ 1 I L/& (2s) .(3.11) 



This integral is identical with that of the bulk correlation function (|2.15 ), 
except that the argument x/a in ( 2.15|) is replaced here by L/a. Therefore 
the large-L behavior of ("|3. 1 1|) is analogous to the large- la; I behavior of fl2.16|), 



An'i-d 

D (x~\ L, a) = ^-^nLraf^ 2 [smH~a/^)f- 3)/2 e "V6 [i + 0(L _1 

(3.12) 



for L ^> £x where £i is the exponential bulk correlation length in the direction 
of one of the cubic axes as determined by (|2.18 ). The exact parallelism 



between the large- distance behavior of the correlation function C(x) and 
the large-L behavior of the susceptibility x{^i L) is the central result of this 
paper. On physical grounds it is quite plausible that for L ^> £ there exists 
a sensitivity of finite-size effects to the length £i governing the large-distance 
decay of C(x) rather than to an averaged length as represented by the second 
moment £ of the correlation function. 



The result ( |3.12| ) can be extended to a d- dimensional system with partially 



finite geometry that is confined in d dimensions and infinite in d — d di- 
mensions. In this case the prefactor d in ( 3.12Q is replaced by d. For 



d — 1 (film geometry) we find agreement with the result of Barber and 
Fisher [[31], |32|. The authors did not recognize, however, that their quantity 
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Ta{T) = 2 arsinh {^q 2 /2j is identical with the inverse of the exponential 
bulk correlation length £i/a. 

The result ( |3 . 3| ) together with ( |3.12j ) is still valid for arbitrary £i/a even 



well above T c . Using the known expression for the bulk susceptibility Xb [10 



we obtain from fl3.3|) and (|3.12|) the relative deviation from the bulk critical 



behavior for L > (i > o as 



A X ^ =9{L/Ci) (3-13) 

Xb 



with the universal function in the large-n limit for 2 < d < 4 



g(L/Ci) = 2 d vr 1 / 2 [r ((4 - d)/<l)]- x (2^/Lf' 1)/2 e . (3.14) 

This result agrees with and goes beyond our previous result in (132) - (134) 
or (135) of Ref. |T0[ which was expressed in terms of £ rather than £i (com- 
pare also ( |1.4p and ( |1.1(J| ) of the present paper). Previously we did not yet 
recognize the physical origin of the non-scaling contribution R(L/£,a/£) in 
(134) of Ref. fTOfl . Now we see that £ e // as defined in the paragraph after 



(107) of Ref. |TIJ turns out to be identical with £i. This is parallel to the 
bulk order-parameter correlation function of Section 3. Thus our previous 
interpretation in terms of a violation of finite-size scaling below four dimen- 
sions in the region L ^> £ was incomplete for the lattice system (but not 



for the continuum system with a sharp cutoff |T0| , 29, 33|). In the critical 
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region £ 3> L, on the other hand, the natural reference length remains to 
be £ and not £1. In this region the function has a power- law form (not 
exponential) in which the non-universal part of the difference between 
and L/£ can be considered as a negligible correction. 

We conclude that nonuniversal and non-negligible lattice effects do exist in 
the region L>( but they can be absorbed in the finite-size scaling argument 
by employing the exponential bulk correlation length. This remedies the ap- 
parent violation of finite-size scaling found previously below four dimensions 



|flO| and simplifies the physical picture of critical behavior in confined lattice 
systems with periodic boundary conditions. Nevertheless we maintain that 
a scaling description of the entire region < L/£ < oo requires to embody 
in the scaling function a kind of crossover in the finite-size scaling variable 
from x = L/£ for < x < 0(1) to x — for 0(1) < x < oo. 

So far our conclusions have only been shown to be correct in the spherical 
(large-n) limit for 2 < d < 4 (and in one-loop order, see Section 4, see also 
Section 6 for the d — 1 Ising model). We note that Eq. fl3.14| ) has a finite 



limit also for d — > 2 at fixed £i. Further work is needed to prove whether 
finite-size scaling for lattice systems with periodic boundary conditions be- 
low four dimensions is indeed an asymptotically exact property for finite n 
beyond one-loop order. General renormalization-group arguments |3], [| are 



not sufficient for such a proof, as shown in Sect. 3 of Ref. |T0|| . 



For a corresponding analysis of the large-L behavior above T c within the 
continuum ip 4 theory we refer to Refs. [T(| |29|, |33| . In this case an exponential 



size dependence of Ax is found only for a smooth cutoff whereas a sharp cutoff 

23 



implies a nonuniversal non-exponential L dependence of Ax |TD 
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4 Perturbative treatment of finite-size effects 
for d < 4 

In this Section we present two different perturbative treatments of the finite- 
size effects of the lattice model ( |2.1[ ) for finite n. We shall focus our interest 
on A% in the region L ^> £ above T c where lattice effects are expected to 
be non-negligible according to the exact results of the preceding Section. 
In particular we show that only the first version of the perturbative treat- 
ment (in Subsection 4.1) correctly predicts the exponential size-dependence 
of Ax oc e~ L/ ^. 



4.1 Ordinary perturbation theory 

First we use ordinary perturbation theory with respect to u Q without sep- 
arating the lowest (k = 0) mode of <£>(x). In one-loop order above T c the 



inverse (bare) susceptibility of the lattice model (|2.1|) in a cubic geometry 
with periodic boundary conditions is given by pi 



1 + 4(n + 2) u J 2 eD (£, L, a) + O(u 2 ) 



(4.1) 



with 
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D = L- d j2 (r 2 + ki Jo 



1 'c' + k/JoY 1 (4.2) 



where £ is the second-moment correlation length. The function 



D(£,L,a) = J D(J r,L,a 



(4.3) 



can be represented in terms of Bessel functions according to ( 3.10Q and ( p. 11 
Eqs. (|4.1|) - (|4.3|) are valid for general d > 2. Because of the k = term in the 



sum of (|4.2|) the perturbative expression (|4.1|) is not applicable to the region 
£ ^> L. In this region a separation of the lowest mode from the higher modes 
is necessary (see Sect. 4.2). But here we are interested in the region L ^> £ 
where the function D is well behaved according to ( |3.12| ). Applying the RG 
procedure of Ref. to the bare expression ( |4.1|) and using the asymptotic 
form (|3.12|) leads to the scaling result for d < 4 and for L ^> ^ ^> a 



A X = g(L/^) = 4(n + 2)u* d (2vrL/ei) (1 ^ )/2 e ~ L/ ^ + 0(u* 2 ) (4.4) 



where u* is the fixed point value of the renormalized coupling [Rll and where 
£i is the exponential bulk correlation length given by ( |2.18[ ), up to two-loop 
corrections. Eq. (|4.4|) has the same form as ( |3.14|) . It also agrees with 



and goes beyond our previous perturbative result (106) of Ref. [ITU]. Now 
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we see that it is £1 = £/[l — (a/£) 2 /24 + ...] rather than £ that should be 
employed in the scaling representation, similar to the case n — > oo discussed 



in the preceding Section. Thus the interpretation of Ref. |10j in terms of a 
violation of finite-size scaling was incomplete since the lattice constant a can 
be absorbed in ^ in a natural way. The conclusions drawn in the preceding 
Section after ( 3.14j) regarding the validity of finite-size scaling in terms of 
£1 apply also to finite n, at least in one- loop order. (Beyond one- loop order 
we expect that the exponential part of (|4.4| ) contains the bulk correlation 
length £1 whose universal amplitude ratio lim 4 ^o+ is slightly larger than 



1 |L4], |l6fl.) Clearly these conclusions can be extended to a <i-dimensional 
system with partially finite geometry that is confined in d dimensions and 
is infinite in d — d dimensions. In this case the result ( |4.4| ) remains valid 
except that the prefactor d should be replaced by d fllUfl . The result of this 
Subsection will be extended to the case of the continuum ip A theory in a 



separate paper 29 



4.2 Separation of the lowest mode 

In the following we discuss the result for Ax for L ^> £ if the standard finite- 
size perturbation theory |TS|, |20], |2T], [32| is used. The details of the calculation 
are given in Appendix C. In this approach the lowest mode is separated 
and treated exactly whereas the higher modes are treated perturbatively. 
Accordingly we decompose 
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(4.5) 



and H = H + H' with 



//„(*) = L d [ ir $ 2 + n $ 4 j . (1.6) 



The susceptibility ( |3.1j ) is expressed as 



x =-/^) = — /^$$ 2 P($) (4.7) 
n n J 



where 



P($) = exp [-# e// ($)] / J d n $ exp [-# e// ($)] 



(4.8) 



is the order-parameter distribution function with the effective Hamiltonian 



if e //($) = # ($) + r ($) . (4.9) 

The present approach consists of a perturbative expansion of T ($) in the 
exponent of P{$) and not of an expansion of x itself. It turns out that this 
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approach does not capture the correct (exponential) size dependence of A% 
for L 3> £ but instead yields 



Ax oc L~ d (4.10) 

in any finite order of perturbation theory (see Appendix C). 

The failure of this approach is due to the fact that the separation of the zero 
mode |TP|, |2I], [H], is inadequate in the region L ^> £. In this region the 



zero mode does not have a dangerous character and all modes including the 
k = mode should be treated in the same way. This argument is valid for 
general d > 2 including d > 4 (Sect. 5). The amplitude A(uq) of the spurious 
power law Ax oc A(uo)L~ d is only partially cancelled order by order in a 
perturbative treatment of the higher modes but A(uq) remains nonzero at 
any finite order of perturbation theory. A complete cancellation of A{uq) 
is achieved only in an exact treatment of the k ^ modes as can be seen 
from the exact solution for n — ► oo (see Eqs. (12) and (21) of Ref. pfl). 
These considerations are insensitive to the lattice spacing and remain valid 
also within the continuum <y? 4 theory |29 . 



We conclude that the perturbative calculation of finite-size effects above and 
at T c requires two different approaches depending on whether < L/£ < 
0(1) or 0(1) < L/£ < oo. In the former case the separation of the low- 
est mode is appropriate. In the latter case which includes the approach to 
the bulk limit at fixed T > T c one should employ the ordinary perturbation 
approach of Subsection 4.1 where all modes are treated perturbatively. To 
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combine the results of both approaches requires some matching in an inter- 
mediate range of L/£. The good agreement of our previous finite-size calcu- 



lations |21j with highly accurate Monte-Carlo data for the three-dimensional 
Ising model [|35| was restricted to the region < L/£ < 0(1) whereas the 
region 0(1) < L/£ < oo was not investigated. The exponential size depen- 



dence in the latter region is not correctly included in the results of Ref. [21 



The same criticism applies to other finite-size calculations in the literature 
which are based on the separation of the lowest mode. Although these expo- 
nential effects are small they are detectable and clearly distinguishable from 
power-law terms as has been demonstrated very recently by Monte-Carlo 
simulations for the magnetization of the two- and three-dimensional Ising 
model p6 |. 
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5 Finite-size effects for d > 4 



In this Section we extend our study of finite size effects to d > 4 within 
the lattice model fl2.1| ) for cubic geometry and periodic boundary conditions. 
We focus our interest on the approach of the susceptibility x to the bulk 
susceptibility Xb above T c . To provide a correct description of this approach is 
a basic task of finite-size theory. This corresponds to the region L ^> £ where 
the exponential correlation length £1 is expected to become an important 
length scale according to the results of the preceding Sections. We shall show 



l/(d-4) 




that, in addition to £1, the second reference length |22|, ^7], |3~6[1 lp ~ u 
associated with the higher (inhomogeneous) modes remains relevant for the 
large-L behavior of x an d that a single-variable (lowest-mode) finite-size 



scaling description 19, 24 of x f aus f° r L ^ £• 



5.1 Exact results for n — > oo 

For x b 3> a 2 the inverse bulk susceptibility above T c for 4 < d < 6 and 
n — > oo at fixed uqu is determined by 



X^ 1 = r„ - r 0c - 4 Mo n x 6 _1 J?\l + (d - A)' 1 fa 1 a 2 ) 



\-l ^.-l%2\(d-4)/2 



(5.1) 



as follows from ( |3.3| ) for L — > 00. From ( |5.1|) and ( |3.3| ) we then obtain the 
leading relative deviation of x from Xb for Ax C 1 as 
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Ax 



Xb-X 
Xb 



If XbJoD(xt\L,~a) + 0[(A X ) 2 } 



(5.2) 



with the reference length [^7 



-I V(d-4) 



Au Q n 



J 2 ( 1 + 4u n f J k : 

k 



(5.3) 



The bulk susceptibility can be expressed in terms of the second-moment cor- 
relation length ^ as Xb = Jq 1 £ 2 - Using the large- L behavior of D(xb ■> L, d) 
according to ( |3.12| ) we obtain from ( |5.2j ) for L ^ £ 



,(l-d)/2/ 



(L// ) 4 - rf (L/e)- 2 [(L/a) sinh(a/Ci)] 



(d-3)/2 



(5.4) 



where now the exponential bulk correlation length £i, (|2.18| ), governs the 
exponential size dependence, similar to the case d < 4. In the asymptotic 
region (>awe may replace £ by £i in the non-exponential part of ( |5.4|) . 
This yields the two-variable finite-size scaling form 



A X = g(L/£ 1 ,(L/l ) 4 - d ) (5.5) 
with the exact scaling function for L ^> £ 3> a 
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g(x, y)= d (2vr)( 1 -^/ 2 y x^ 12 e~ x . (5.6) 

Unlike the corresponding scaling function g(x) in (|3.14j) for d < 4 we see that 
here we need two scaling variables x = L/^\ and y = (L/lo) 4 ~ d . In the present 
context where the lowest mode plays no particular role, the second variable y 
is associated with the higher modes and has nothing to do with the dangerous 
character of Uq. The present result ( |5.5| ) and ( |5.6| ) complements our previous 
two- variable finite-size scaling function (138) - (142) for the lattice model in 



Ref. pTj where £ instead of £i was employed |37|]. A complete description of 
the scaling form of x i n the entire (asymptotic) L^ 1 — plane requires to 
incorporate in g(x,y) a kind of crossover from the variables (L/£,y) for the 
region < L/£ < 0(1) to (L/^,y) for the region 0(1) < L/& < oo. 

We recall that an alternative choice of the scaling variables (L/£,y) is (w,y) 



where [22 



w = (L/0 2 y- 1/2 = t(L/if 2 , I = / (W/o) 4/d • (5.7) 



Correspondingly the susceptibility can be represented as [22 



X = L d ^P(w,y). (5.8) 
Instead of (w,y) an equivalent choice is (w 2 ^ d ,y) where 
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w 2 ' d = L/£ 7 



(5.9) 



contains Binder's "thermodynamic length" It |38|] which is related to £o and 
Zo as 



i T = if-^ d f/d = if~^ d £ /d r 2 / d . (5.10) 

This length scale, together with Iq, plays an important role in the region 
< L/£t ^ 0(1) where the dangerous character of uq is important (cor- 
responding to the region between the curved lines in Fig. 1 of Ref. p2| ) 
but £t loses its significance outside this region. In particular in the region 
0(1) < L/£ < oo the correlation length £i (and £) and the reference length Iq 
associated with the higher modes [p| govern the finite-size effects, as demon- 
strated by (|5\4l ) - ( p^6|) . Thus, not £t alone but £i and / are indispensable 
for a complete description of the finite-size effects in the entire asymptotic 
(large L, small t) region for d > 4. In the following Subsection we show that 
ignoring the lengths £ or £x and Iq implies an incorrect large- L dependence 
of x at any fixed T > T c . 



5.2 Lowest-mode approximation for n — > oo 

Neglecting the k ^ contributions to \ leads to 
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Xo = 2 



d\V2 



-i 



(5.11) 



for n — > oo at fixed u n. In this approximation we have X° = r o 1 = ( a o^) 1 - 



This yields 



Axo^^^ = 1-2 



1 + (1 + 16u Q n r^ 2 L' d ) 1/2 (5.12) 



In the present approximation the lengths £o and Iq are reduced to 



£o = (J /ao) 1/2 , k = (4u ri J 2 ) 



2 \l/(d-4) 



(5.13) 



Thus Axq can be expressed in terms of the thermodynamic length Q5.1U ), 



as 



A Xo = 1-2 1+ 1 + 4(VI) 



(5.14) 



(L/£ T )- d + 0(L- M ) 



(5.15) 



for L 3> ^t- Comparison of ( J5.14Q with ( 5T1 ) - ( |5.6| ) shows that the lowest- 
mode approximation [19] fails both with regard to the L dependence of Ax 
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as well as with regard to the temperature dependence of the reference length 
scale It ^ £1 in the scaling variable. 



The same criticism applies to the phenomenological single-variable scaling 
form A% = /(L/£y) proposed by Binder et al. G4l]. The recent statement 



39fl that the single-variable scaling form is presumably true asymptotically 
for L — > oo is correct for T = T c (more precisely, for L —>■ oo at fixed finite 
w). Furthermore, the lowest-mode result Xo{t, L) = L d ^ 2 P(w,0) correctly 
contains the limit Xb(t) = Xo(^ 00 ) of Xo(t,Li) for L — > 00 at fixed L/£ (see 
Eq. (104) of PH and Eq. (102) of [0) but does not correctly describe the 
size dependence in approaching this limit Xbft) ( see -^l- (104) of [0). A 
corresponding statement is also true with regard to the magnetization below 



T c as confirmed by Monte Carlo simulations for the d — 5 Ising model |40 
As pointed out in pi, the property 



lim X (t,L)/ Xb (t) = 1 (5.16) 

at arbitrary fixed L/£ is a nontrivial feature that is valid only for d > 4 
and that is correctly contained in the lowest-mode approximation and in the 
Binder et al. scaling form whereas for d < 4 the same limit yields the function 
/(L/C) ^ 1 for LH < 00. 

For any fixed T 7^ T c , however, the leading size dependence oc L~ d predicted 



by the single- variable scaling form |24] is incorrect. The origin for this defect 



are the missing higher modes. At fixed w for large L, these modes only 
cause (slowly decaying) corrections ~ 0(y 1 ^ 2 ) to the leading size dependence 
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L d / 2 P(w,0) of the lowest-mode approximation [22"]. For fixed T ^ T c and 
large L corresponding to w ^> 1 and y <C 1, however, the higher modes and 
the lowest mode must be treated in the same way, as shown in Sect. 4.2, and 
the effects of the higher modes become increasingly dominant with increasing 
L/£ and can no longer be considered only as corrections. 

More specifically, the structure of the scaling function (|5.8|) can be written 

as 



L d ' 2 P{w, y) = L d ' 2 P{w, 0) + A{t, L) (5.17) 

where A describes the size effect of the higher modes. The size effect Ao of 
the zero mode is contained in 



L d ' 2 P(w,0)= X b(t)+A (t,L) (5.18) 

where Xb(t) is the bulk susceptibility. The crucial point now is that for 
sufficiently large L ^> £ the structure of A becomes |33| 



A = -A + O(e- L/?1 ), (5.19) 



i.e., the zero-mode size dependence of A oc L d is exactly cancelled by the 
higher- mode size dependence of A. Thus, A is not small compared to A 
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and it is inadequate to refer to the size effects of A only as "corrections to 
the lowest-mode result" |39|, [4l| . 



5.3 Perturbative treatment for finite n 

For finite n a perturbative treatment of the finite-size effects becomes nec- 
essary. Our arguments (in Sect. 4) for the necessity of two different per- 
turbative approaches remain valid also for d > 4. A one-loop perturbation 
calculation on the basis of a separation of the lowest mode was presented 
recently |22 for the case n = 1. The results of this calculation are applicable 



to the region < L/£ < 0(1) but the quality of this approach deteriorates 
with increasing L/£ in the region 0(1) < L/£ < oo. For the latter region the 
following ordinary perturbation calculation with respect to u is appropriate. 

The bare perturbative one-loop expression ( [4.1| ) remains valid also for d > 4. 
This leads to 



A X = 4(n + 2)u J 2 e D(£, L, d) + O(u 2 ) (5.20) 



where D is given by ([4. 3|) and ( |3.10 ). Here we interpret the prefactor 4(n + 



2)u Jq 2 as Zq -4 + 0(ul) as indicated by the result ( |5.2| ) for x i n the large-n 
limit. Using ( |4.3|) and ( p,12| ) for L ^> ^ we arrive at the same expression for 
Ax as given already in ( |5.4|) where now the reference length Iq for finite n is 
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Au (n + 2) 



l/(d-4) 



I 2 



1 +Au (n + 8) J J k 



(5.21) 



The coefficient 4(n + 8) in the denominator is inferred from the form ( |5.3| ) 
in the large-n limit and from the previous result in (30) of Ref. |22j for 
n = 1. (Our present definition of /o differs from that of Ref. |22]| by the 
factor [4(n + 2)] 1 ^ d ~ 4 ^). In the asymptotic region L > ( > a we arrive at 
a one-loop scaling form of Ax for finite n which is identical with ( |5.5|) and 
(|5.6|). The n dependence enters only the expressions for l Q and 



£o = (Jo/«o) 1/2 



1 + 4(n + 2)u / J, 



1/2 



(5.22) 



For n = 1 this result for A% complements our previous two-variable finite- 
size scaling function (97) and (99) in Ref. which did not yet incorporate 
the exponential size dependence oc e~ L ^ x for L 3> 



Similar to the case n — > oo, the lowest-mode approximation |19| and the 



phenomenological single- variable scaling form |24| lead to the incorrect power 
law for finite n 



+ 0(L- 2d ) = (L/i T )- d + 0(L-' 2d ) 



-2d\ 



(5.23) 
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with i T given by ( |5.1(J| ). Our comments after ( |5.14| ) apply also to ( p.23| ). 



In summary, while the two-variable scaling form correctly embodies the in- 
dependent approach oc e~ L ^ x to the bulk limit at fixed T > T c this crucial 
information is lost in the single- variable scaling form and in the lowest-mode 
approximation. Therefore the reduction of the two scaling variables to a 
single scaling variable is not justified. 

Very recently our prediction of the exponential size dependence ( |5.6| ) has 
been confirmed by Monte Carlo simulations for the five-dimensional Ising 



model 25, 261. 



An analysis of finite-size effects for d > 4 within the continuum (p 4 theory 
will be given in a separate paper [^{J where it is shown that the finite-size 
effects depend significantly on the cutoff procedure. For a smooth cutoff 
the results for the Binder cumulant at T c and for the two- variable finite-size 
scaling function of \ are found to be different from those found previously 
22], |27|, [33], for a sharp cutoff, see also the note added in Ref. |10| . 
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6 One-dimensional Ising model 

In this Section we illustrate in an elementary way the close connection be- 
tween the exponential bulk correlation length £1 and the finite-size scal- 
ing structure for the example of the exactly solvable one-dimensional Ising 
model. Although the critical temperature T c = vanishes, this model has 
well defined correlation lengths £ and (i / ( for T > 0, which diverge for 
T -f T c = 0. 

First we consider spins Sj = ±1 on sites Xj of an infinite chain with a lattice 
spacing a. The Hamiltonian reads 



H — J s j s j+i 



(6.1) 



j=-oo 



The correlation function is well known, see e.g. Ref. The exact result 

has an exponential form for arbitrary distances \xi — Xj\, 



< SiSj >= exp(-|xi - a?j|/£i) (6.2) 
with the exponential correlation length 



6 = a [ln(A + /A_)] 



-i 



(6.3) 
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where A + and A_ are the eigenvalues of the transfer matrix [53] with A + > A_. 
Obviously (|6.2| ) has a scaling form in terms of £1. 



In order to calculate the correlation length £ as defined in (|1.2| ) we consider 
the Fourier transform 



oo 

G(k) = < s Sj > exp(— i k j a) (6.4) 

j=-oo 

where so denotes the spin on a fixed site Xq. Using ( |6.2| ) and ( |6.3|) we obtain 



G(k) 



1 ~ (A^/A + ) 2 
1 + (A_/A+) 2 - 2(A_/A+) cos A; 



(6.5) 



This leads to the exact result G(0) = exp(2/3e) and 



e 2 = G(0)A G{k) 



-1 



a 2 



4 [sinh(a/2£i)]' 



(6.6) 



where £i is given by ( |6.3|) . This relation between £ and £x is identical with 
( |2.17| ) or ( |2.18| ) which was derived for the (p A model in the large-n limit in 
Sect. 3 and in one-loop order in Sect. 4. In particular we again have — > 1 
for T T r = 0. 



Now we consider the finite-size effect on the susceptibility xl of a finite one- 
dimensional Ising chain which consists of iV spins and which has a length 
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L = Na. We assume periodic boundary conditions. The partition function 
43| Zn = + X N which is valid also at finite magnetic field h. The 



is 



second derivative with respect to h leads to the exact finite-size scaling form 
of the relative deviation from the bulk susceptibility at h — 



A Xoo Xl 2 e Z^ 1 . . 



Eq. ( |6.7|) is valid for arbitrary where is identical with the exponential 
bulk correlation length (|6.3|). 

The crucial point is that Z//£i and not L/^ is the natural finite-size scaling 
variable. If Ax were expressed in terms of £ then an apparent violation of 
finite-size scaling would arise from the a dependent difference between £1 and 



A X = 2 e~ m exp (La 2 /24£ 3 ) (6.8) 



for L ^> £ ^> a, in the same way as found previously JTOf for the y? 4 model. 
Similarly, the bulk scaling form for the correlation function would be violated 
if the result ( |6.2| ) would be expressed in terms of £. Thus the exact results 
(6.2), (|6.6| ) and ( |6.7|) demonstrate in an elementary way the significant dif- 



ference between £ and (i as well as the close connection between bulk and 
finite-size scaling. 

It would be interesting to extend this analysis to the exact results for the 
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two-dimensional Ising model [13, 35] and to compare these results with the 
exponential size dependence found in recent Monte Carlo data in Fig. 2c of 
Ref. ||. 



Note added 

The distinction between £ and £i is significant also for resolving discrepan- 
cies in the interpretation of Monte Carlo simulation results of percolation 



phenomena |f46|. 
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Appendix A : Anisotropy of the exponential 
correlation length 



We start from (|2.14j ) for the case where x = (x%, X2, 0, ...), 



a 2 ~ d 



G(x) = \j~ I ds e~^ s e- 2ds [h{2s)] d ~ 2 I Ui {2s)L 2 {2s) (A.1) 



with Pi = \xi\/a. For large v and large s = uz/2 we have [17 



I v (vz) ~ (27rz/)- 1 / 2 g- 1 / 2 exp (z/ {g + In [z(l + q)' 1 ] }) (A.2) 



with q = (1 + .2 2 ) 1 / 2 . Furthermore we use the large-s behavior [^ 



/ (2s) =e 2s (4ns)- 1/2 [l + O^ 1 )] . 



(A.3) 



For sufficiently large a saddle-point approach suffices to perform the inte- 
gration over s and to determine the exponential large- \xi\ behavior of ( |A.l ). 
Introducing the angle 9 according to v x = rcosO, v 2 = r sin 8 we obtain the 
exponential part of G(x) as 



G(x) ~ exp {-\x\M8)} 



(A.4) 



45 



where denotes the anisotropic exponential correlation length. The an- 
gular dependence is given by 



= (cos 6) In [u 1/2 cos0 + (1 + wcos 2 £) 1/2 ] 

+ (sinS) In [u 1/2 sin^ + (1 + u sin 2 6) 1/2 ] (A.5) 



where 



u 



b(b 2 -A) b+ (b 2 sin 2 29 + 4 cos 2 29) 



1/2 



-i 



(A.6) 



6 = 2 + (S/072 



(A.7) 



For the case 9 = this yields 



6(0) 



arcsinh 



2s c 



-i -l 



(A. 



corresponding to (|2.18|). For the case 9 = 7r/4 the result is 



&0r/4) 



2 3/2 



arcsinh 



2 3/2^ 



(A.9) 



46 



Eqs. (|Op and ( |OD agree with (4.12) and (4.14) of |T3| for d = 2. For 
(»awe obtain from ([OP - (|OD 



6W e 



l-±(lW26>)(|) +0(aVe 4 ) 



(A.IO) 



For = and 9 = vr/4, (jA~TTD disagrees with (4.13) and (4.15) of Ref. 
|13| . Asymptotically (£ — > oo) we obtain from (|A.10j ) — > 1 for all 6 
(in the spherical limit and in one-loop order). Nevertheless, because of the 
exponential form of G(x), the nonasymptotic ^-dependence remains non- 
negligible in Cr(x) for |x| > 24£ 3 /a 2 even arbitrarily close to T c , see the 
reasoning in the context of ( 2.22| ) and ( p. 23 ). 



To derive the large-|x| behavior (|2.16|) for the case x = (x, 0, ...) we use (|A.2|) 
and ( |A.3| ). Expanding around the maximum of the exponential part of the 
integrand of C(x), (|2.15|) , and performing the integration over s leads to 
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Appendix B : Bulk correlation function in one- 
loop order for d < 4 



In terms of the second-moment correlation length Q1.2Q the bare bulk two- 
point vertex function at finite k above T c is given in one-loop order for d > 2 

by m 



(k, £, u , a, d) = J k + J C 2 + 0{ul) 



(B.l) 



corresponding to the integrand of G(x) in the form of ( [2.11 ) and Q27HX ). Em- 
ploying the renormalization procedure at finite lattice constant a as described 
in Section 2.2 of Ref. []IU] we obtain for d < 4 



r( 2 )(k) = Z~ l J k + J r 2 + O (u(£) 2 ) exp / C v (0^- • 



(B-2) 



For the application to T > T c we choose the flow parameter as I = £o/£ where 
£o is the asymptotic amplitude of £ above T c . For the case x = (x, 0, 0...) and 
in the limit \x\ ^> a the correlation function is in one-loop order for d < 4 



G(x) 



[r( 2 )(k 



ikx 



C(x,6,a)exp y C P (^j (B.3) 
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with C(x, £i,a) given by ( |2. 16| )-( PTT5| ) , apart from corrections of 0(u(£) 2 ), 
where u(£) is the effective four-point coupling jlOfl . In the asymptotic region 
(1 > a we obtain from ( |B.3| ) 



G(x) = Z v [A^] 1 r" Co C(x, 6, a) [1 + O (u* 2 )] (B.4) 

with 



A( 2 ) = exp{^°[C,(0-aO)]f } 



(B.5) 



and the critical exponent 77 = — (^,(0). In ( |B.4|) the asymptotic (\x\ ^> £1) 
form of C (x, £1, a) is 



~ n 2-d / ~ \ (d-i)/2 / ~ X (d-3)/2 

c <^-)-V(4?) ! e " w *- (R6) 



Asymptotically we may replace ^ in ( |B.4j ) by £1 [1 + (£ 2 )]- Therefore 
G(x) can be rewritten in the asymptotic scaling form for d < 4 



G(x) = (a/|x|) rf - 2+ " $(|x|/6) (B.7) 



with the scaling function for |x|/£i 1 
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*(M/6) = A u ^-m (^) 2( ^^(-H/eo, 

apart from O (w* 2 ) corrections. The amplitude A is for d < 4 []T0| 

I = z v («, a/e , d) Ceo/a)' [aW]" 1 . 
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Appendix C : Separation of the lowest mode 



In the following we show that the separation of the lowest mode implies 
A% oc L~ d for L>( and for general d > 2 in any finite order of perturbation 
theory. In leading order of the k ^ modes the effective Hamiltonian 



becomes [22 



H eff ($) 



eff 



u, 



eff 



$1 



(C.l) 



.eff 



r - r 0c + 4(n + 2)u 



L- d ( r o - r oc + Jk 



4" 1 



, (C.2) 



u 



eff 



Wo - 4(n + 8) ul L d ^ ( r o ~ r 0c + X 



(C.3) 



Since we work here at finite lattice spacing we have incorporated the finite 
shift r 0c = — 4(n + 2)u J J^ 1 + O (uq) of the parameter r already at one- 

k 

loop order. For simplicity we have not included here the zero-mode average 
Mq defined previously |2T|, |22| since it is negligible in the region T > T c . It 
is convenient to rewrite Tq and in terms of the second-moment bulk 



correlation length £. Using the bare one-loop relation [1C] 
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r -r Qc = J £ 2 < 1 + 4(n + 2)u 



Jk [ Jk + Jo £ 



-2 



+ O (ul) 



we obtain 



(C.4) 



J r 2 - 4(n + 2) n J^ 1 A x (£~ 2 ) 



(C.5) 



u 



eff 



U 



4(n + 8)u 2 J 2 / (c 2 + jjJoY +4(n + 8)u 2 J 2 A 2 (C 2 ) , 



(C.6) 



A m (r 2 ) = / ( r 2 + A/ j, 



(C.7) 



Note that, because of the separation of the zero-mode, the sums in (|C.6|) and 
( |C.7|) do not contain the k = part that was still contained in the function D, 
(gg ). Therefore we obtain the relation A x (£~ 2 ) = L- d £ 2 - L, a). The 
important consequence is that the power- law term oc L~ d in Ai(£ -2 ) becomes 
dominant for large L/£ compared to the exponential behavior oc e _L ^ 1 of D 
according to ( [4.3| ) and ( |3.12j ) . Similarly we have the non-exponential behavior 
A2 oc L~ d for large 
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In the present approximation the susceptibility is 



L d n 1/2 



n \u 



X=-[- m ) A (Y^) , (C.8) 



Y°" = LW r tf f (u ff )' 1/2 , (C.9) 



/ ds s n+1 exp (— \ y s 2 — s 4 ) 

Mv) = ^ • (c- 10 ) 

J ds s n_1 exp (— | y s 2 — s 4 ) 
o 

These expressions are valid for general d > 2 in the sense of bare perturbation 
theory. If we take the limit 5^0 and apply the field-theoretic RG approach 
to (|08D - (ICTTOl) we reproduce the results of Ref. J2TJ for 2 < d < 4 which 
are well applicable to the critical region 0<L/£<O(l). In this region 
lattice effects are negligible. 

In the region L/£ ^> 1 corresponding to Y e ^ ^> 1 the function d 2 has the 
asymptotic expansion 



Mv) = n y- 1 [1 - 4(n + 2)y- 2 + O (y^)} . (C.ll) 

Within the lowest-mode approximation, where = Jo£~ 2 an d u e J^ = u , 
the expansion ( |C.11| ) implies 
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X = Jo 1 e [1 - 4(n + 2) n J - 2 ?/L d + O (L~ 2d )] , (C.12) 



A X = 4(n + 2)u J 2 e 4 L~ d + (L~ 2d ) . (C.13) 

We see that the lowest-mode approximation yields the incorrect (non-expo- 
nential) size dependence oc L~ d for large L. In order to see whether this 
defect is remedied in higher order we proceed by including the next terms of 
( |C.5|) and ( C.6 ). Then it turns out that the 0(uq) contribution oc Ax oc L~ d 



in ( |U.5| ) cancels the O(u ) term in ( |C.12j ) and ( |U.13| ). Therefore A% appears 
to become proportional to according to 



A X = -cu 2 Jo 4 Z 8 - d L- d + 0(L- 2d ) , (C.14) 
c = 8(n + 2)(n + 8) A d (d - 2) (4 - d)~ l , (C.15) 
A d = T(3 - d/2) 2 2 - d n- d/2 (d - 2)- 1 , (C.16) 

for 2 < d < 4. But in this order the leading size dependence Ax oc L~ d in 
( |C.14|) is still non-exponential. We anticipate that in the next order of this 
approach the Uq term of ( |C.14| ) will also be cancelled, thus A\ oc ulL~ d , etc. 
We conclude that the present approach does not capture the correct (expo- 
nential) size dependence of Ax for L ^> £ in any finite order of perturbation 
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theory. An application of the renormalization group for d < 4 would not 
remedy this defect for L ^> £. In the present context the renormalization of 
(|C14 ) would only change Jq 2 uq to u* £ d ~ i and this would only change the 
critical temperature dependence of ( |C114j) to A% = —cu* 2 (L/£ t )~ d for d < 4 
in the region L > ( > o. 
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